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Abstract. We consider the problem of efficient integration of an n-variate polyno- 
mial with respect to the Gaussian measure in R n and related problems of complex 
integration and optimization of a polynomial on the unit sphere. We identify a class 
of n-variate polynomials / for which the integral of any positive integer power f p 
over the whole space is well-approximated by a properly scaled integral over a ran- 
dom subspace of dimension O(logn). Consequently, the maximum of / on the unit 
sphere is well-approximated by a properly scaled maximum on the unit sphere in a 
random subspace of dimension O(logn). We discuss connections with problems of 
combinatorial counting and applications to efficient approximation of a hafnian of a 
positive matrix. 



We consider the problem of efficient integration of multivariate polynomials with 
respect to the Gaussian measure in K n . 

Let us assume that the real n-variate homogeneous polynomial / of degree m is 
given to us by some "black box", which inputs an n- vector x = (£1, . . . , £ n ) an d 
outputs the value of f(x). We want to compute or estimate the integral 



(27r)-"/ 2 e-l |a; l |2 / 2 , where ||x|| = £ + ■ ■ ■ + £ for x = (£1, . . . , £ n ). 
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1. Introduction 




where \x n is the standard Gaussian measure with the density 
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Typeset by Aj^S-T^i. 



If to is odd then the integral is 0, so the interesting case is that of an even degree 

TO. 

An equivalent problem is to integrate / over the unit sphere S n_1 C M n . As- 
suming that to = 2k is even, we have 

/ f( x ) dx = okrt^/o! lA f f 
J §n -i 2*T(n/2 + k) J Rn 

where dx is the rotation invariant Haar probability measure on § n_1 . This and 
related formulas for integrals of polynomials over the unit sphere and over the 
Gaussian measure on M n can be found, for example, in [B02b]. 

The most straightforward and the most general approach to integration is to 
employ the Monte Carlo method, that is, to sample N random points Xi G 
and approximate the integral by the sample mean: 

1 N 

/ f{x) dx*i — ^2f(Xi). 

Although one can show that for a "typical" polynomial the Monte Carlo method 
works reasonably well, there are simple examples of polynomials where one would 
require to sample exponentially many points to get reasonably close to the integral. 

(1.1) Example. Suppose that f(x) = £f fc for x = . . . , Then 




r(n/2)r(l/2 + fc) 
0Fr(n/2 + k) 



If we choose k ~ n/2 then the integral is of the order of 2~ n for large n. 

On the other hand, if we sample iV random points x\ on the unit sphere S n_1 , 
then with high probability we will have |£i | = 0(y/\n N/n) for the first coordinate £i 
of every sampled point, cf., for example, Section 2 of [MS86]. Thus to approximate 
the integral within a factor c n for some absolute constant c, the number iV of 
samples should be exponentially large in n. 

The reason why the Monte Carlo method doesn't work well on the above example 
is clear: the polynomial f(x) = £f h acquires some large values for an exponentially 
small fraction of x E but those values significantly contribute to the inte- 

gral. In other words, the Monte Carlo method wouldn't work well if the graph of 
the polynomial looks "needle-like". In this paper, we suggest a method tailored 
specifically for such needle-like polynomials. 

The following defines the class of "needle-like" or "focused" polynomials we deal 
with. 
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(1.2) Definitions. Let 

(x,y)=€iTh + ... + € n 7i n for x= (fi,... ,Cn) and y = (7/1, . . . , rj n ) 

be the standard scalar product in R n . 

Let us fix a number < 5 < 1 and a positive integer iV. We say that a ho- 
mogeneous polynomial / : R n — > R of degree m is (5, N)-focused if there exist iV 
non-zero vectors Ci, . . . , cn € R n such that 

• for every pair the cosine of the angle between q and Cj is at least 5; 

• the polynomial / can be written as a non-negative linear combination 

/ iei 

where the sum is taken over subsets I C {1, . . . , N} of cardinality |/| = m and 
ai > 0. 

Our first result is that the value of the integral of a focused polynomial over a 
random lower-dimensional subspace allows one to predict the value of the integral 
over the whole space. 

For a /c-dimensional subspace L C R n , let \ik be the Gaussian measure con- 
centrated on L with the density (27r) -fc / 2 exp{ — ||x|| 2 /2} for x £ L. We pick a 
/c-dimensional subspace at random with respect to the Haar probability measure 
on the Grassmannian Gfc(R n ) and consider the integral 

/ dn k . 

L 

We claim that as long as k ~ log N, the properly scaled integral over L approximates 
the integral over R n within a factor of (1 — e) m / 2 . 

(1.3) Theorem. There exists an absolute constant 7 > with the following prop- 
erty. 

For any 5 > 0, for any positive integer N, for any (S, N)- focused polynomial 
f : R n — ► R of degree m, for any e > 0, and any positive integer 
k > r ye~ 2 d~ 2 ln(iV + 2), the inequality 

(1 - e)™/ 2 / / dfJLk < (-) ' [ f d^ n < (1 - e)~ m ' 2 [ f dfi k 

holds with probability at least 2/3 for a random k- dimensional subspace L C R n . 

Assuming that we can integrate efficiently over lower-dimensional subspaces (see 
Section 1.5 below), we get a randomized approximation algorithm for computing 
the integral of / over R n . Namely, we sample a random /c-dimensional subspace 
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L, compute the integral over L and output the value of that integral multiplied 
by (n//c) m / 2 . To sample L from the uniform distribution on the Grassmannian 
Gfc(lR n ), one can sample k vectors x±, . . . ,Xk independently from the Gaussian 
distribution in R n and let L = spanjxi, ... , xu]- 

One "anti-Monte Carlo" feature of the algorithm is that the estimator is decid- 
edly biased: the expected value of the output is essentially greater (by a factor of 
(n/fc) m / 2 ) than the value we are trying to approximate. This is so because the 
distribution of the integral over a random subspace has a "thick tail": there are 
subspaces which result in large integrals that significantly contribute to the integral 
over the whole space but such subspaces are very rare. 

To increase the probability of obtaining the right approximation, one can use the 
standard approach of sampling several random subspaces and finding the median 
value of the outputs. 

One can observe that if / is (5, iV)-focused then f p is also (5, iV)-focused for any 
positive integer p. This allows us to deduce that the maximum of / over the unit 
sphere is well approximated by the scaled maximum of the restriction of / onto the 
sphere in a lower-dimensional subspace. 

(1.4) Corollary. There exists an absolute constant 7 > with the following prop- 
erty. 

For any 5 > 0, for any positive integer N, for any (5, N) -focused polynomial 
f : R n — ► R of degree m, for any e > 0, and any positive integer 
k > -fe~ 2 5~ 2 ha(N + 2), the inequality 

/u\ m / 2 

(1 - e) m / 2 max fix) < - max fix) < (1 - e)" m/2 max fix) 

holds with probability at least 2/3 for a random k-dimensional subspace L C M n . 

The problem of optimization of a polynomial on the unit sphere has attracted 
some attention recently, see [F04] and [K+04] . Note that by restricting the polyno- 
mial onto a /c-dimensional subspace we effectively reduce the number of variables 
to k in the optimization problem. Using methods of computational algebraic ge- 
ometry allows one to optimize a polynomial over the sphere in time exponential in 
the number of variables. Hence with k = O(logAT), we obtain a quasi-polynomial 
algorithm of rn°^ ogN ^ complexity which approximates the maximum value of the 
polynomial on the sphere within a (1 — e) m / 2 factor. If the degree m of the poly- 
nomial is fixed and iV is bounded by a polynomial in the number n of variables, we 
get a polynomial time approximation algorithm. 

(1.5) On the computational complexity. Let / : IR n — ► 1 be a homogeneous 
polynomial of degree m given by its "black box" which outputs the value of f(x) 
for an input x G M n . Then one can compute the monomial expansion 

/(x) = ^c Q x a where x a = x" 1 . . . x"" for a = («!,..., a n ) 

a 
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in O (( n+ ™ 1 ) 3 j time through the standard procedure of interpolation, cf. also 

[KY91] for the sparse version. If L C W 1 is a /c-dimensional subspace, by choosing 
an orthonormal basis in L, we can identify L with M fc . Then the monomial expansion 
of the restriction f^ can be computed in O (( fc+ ™ -1 ) 3 ) time. If k is fixed, we get a 
polynomial time algorithm. In we choose k = O (log N), the algorithms we obtain 
will be "quasi-polynomial", with the complexity of m°' log N \ 

Once a monomial expansion is obtained, it is easy to integrate polynomials 
since there are explicit formulas to integrate monomials. Given a monomial x Q = 
Xi 1 ■ ■ ■ x% n , the formula is 



/ x a d/j n = I 



n -n/2 Y^ =i 2 ^/2 r (a^i j jf a n a . are even 
otherwise. 



In Section 2, we prove Theorem 1.3 and Corollary 1.4. In Section 3, we con- 
sider some examples and applications, including the problem of approximating the 
hafnian of a positive matrix. In Section 4, we consider the problem of integrating 
polynomials with respect to the complex Gaussian measure in C n . We prove a 
version of Theorem 1.3 in this case and show connections between efficient complex 
integration and certain hard problems of combinatorial enumeration. 

2. Proofs 

One major ingredient of the proof of Theorems 1.3 is the formula for the integral 
a product of linear forms. 

(2.1) Definitions. Let m = 2k be an even positive integer. A perfect matching 
I of the set {1, . . . , m} is an unordered partition of {1, . . . , m) into a union of k 
unordered pairwise disjoint pairs 

I = {{n,ji},{«2,j2},--- ,{«fc,jfc}}- 

Let C = (cij) be an m x m matrix, where m = 2k is an even integer. The hafnian 
haf A of A is defined by the formula 

haf C = c/, 
i 

where the sum is taken over all perfect matchings / of the set {1, . . . , m} and cj is 
the product of all Cy for all pairs {z, j} G /. 

The following result is known as the Wick formula, see, for example, [Zv97] . 
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(2.2) Lemma. Let m be a positive even integer and let ii : M n — > R, i = 1, . . . , m, 
be linear functions. Let C = (cy) be an m x m matrix defined by 



Cij 



= / £i(x)£j(x) dp> n . 



Then 



„ m 

/ Yltiix) dfji n = ha£C. 



If I , is defined by £i(x) = (a^, x) for some G M n then = (a^, a,). 

We also need a version of the Johnson-Lindenstrauss "flattenning" Lemma, see, 
for example, [Ve04] . We present such a version below (with non-optimal constants) , 
taken off Section V.7 of [B02a]. 

(2.3) Lemma. Let x G M. n be a vector and let L C M 71 be a k-dimensional subspace 
chosen at random with respect to the Haar probability measure on the Grassmannian 
Gyt(IR n ). Let x' be the orthogonal projection of x onto L. Then, for any < e < 1, 
the probability that 

(l-6)||x||<y|||x'||<(l-6)- 1 ||x|| 

is at least 1 — 4exp{— e 2 /c/4}. 

The following is a straightforward corollary. We establish it in a slightly larger 
generality than immediately needed, having in mind applications to complex inte- 
gration in Section 4. 

(2.4) Lemma. Let us choose 5 > and e > 0. Suppose that ai, . . . , ajv and 

bi, . . . , 6tv ore vectors from M, n such that the cosine of the angle between every pair 
ai and bj of vectors is at least 5 > 0. 
Let us choose a p > such that 

a-p)- 2 <i + | 

and an integer 

k > minjn, Ap~ 2 In (\2N 2 + 24iV) j . 

Let L C M n be a k-dimensional subspace chosen at random with respect to the Haar 
probability measure on the Grassmannian Gfc(M n ). Let a^bj be the orthogonal 
projection of ai, bj onto L. Then with probability at least 2/3 



(l-e)(a,,6,) < -K,&;.> < (l-e)- 1 ^) 
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for all pairs 

Proof. Scaling, if necessary, we may assume that ||a^|| = ||6j|| = 1 for alH and j, so 
(aj, bj) > S for all We have 

Ktj) JI« + M'-IMI'-INI' and K „ i)= M + W'-W-W\ 

We note that 

(1-P)" 2 <1 + | and (l-p) 2 >l-|. 

Since there are altogether N 2 + 2N vectors a;, bj, and a; + 6j, by Lemma 2.3, for 
a random /c-dimensional subspace L, with probability at least 2/3, we get 



k + y 2 (i - p) 2 < ^ik + &;-n 2 < (i - p)- 2 ik + m 2 



M 2 (1-P) 2 <^K|| 2 <(1-P)- 2 M| 2 and 

2/1 „\2 ^ n ili/||2 ^ /i „\-2||7 1 1 2 



and, similarly, 



! (i-pr <-^;ir<(i-p)- 

for all pairs Since 1 1 cx^ 1 1 = \\bj\\ = 1 and \\a,i + < 2, we get 
hi + bj\\ 2 - ~y < - Wa'i + b'jW 2 < \\a i + b j f + — 

and, similarly, 

.. ..o 5e n.. /ll2 „ n no Se . 

IKH -v-~rKII - °* and 

3 fc 3 

2 -|<^ll 2 <IIM 2 + |- 



Therefore, 



(a,i,bj) -5e< — <a<, 6j> < (a^fy) + 5e 



Since (a;, > 5, the proof follows. □ 

(2.5) Corollary. There exists an absolute constant 7 > with the following prop- 
erty. 

Let 5 > and e > be numbers, let N be a positive integer, and let a\, . . . , ajv 
and b\, . . . , 6/v &e vectors from W 1 such that the cosine of the angle between every 
pair a,i, bj of vectors is at least 5. Let k be a positive integer such that 

k > 7 <rV 2 ln(N + 2) 
7 



and let L C M, n be a k- dimensional subspace chosen at random with respect to the 
Haar probability measure in the Grassmannian GkiW 1 ) . Leta'^b', be the orthogonal 
projections of a^bj onto L. Then, with probability at least 2/3, we have 

(l-e)(a>,b' J }<^(a l7 b J )<(l-e)- 1 (a' l ,b' J ) 

for all pairs a^, bj. 

The proof follows by Lemma 2.4. 

Now we are ready to prove Theorem 1.3. 
Proof of Theorem 1.3. We can write 

f( x ) = ^2aiY[{ci,x), 

i iei 

where the cosine of the angle between every pair of vectors q and cj is at least S, 
I ranges over subsets I C {1, . . . , N} of cardinality m, and oti > 0. For every /, 
let us consider the m x m matrix Ci whose entries c^- are defined by c^- = (q, cj). 
Then, by Lemma 2.2, 

f(x) djji n = a i haf Cj. 



/ 



Let I C R" be a /c-dimensional subspace. Then the restriction f L of / onto L can 
be written as 

h(x) = JJ(c-,x), 
i iei 

where c[ are the orthogonal projections of Cj onto L. Therefore, 



/ f(x) dfi k = J^a/haf C' n 



where the entries of C\ are defined by = (c'^c'j). Since the hafnian of an 
m x m matrix is a no n- negative homogeneous polynomial of degree m/2 in the 
entries of the matrix, the proof follows by Corollary 2.5 where we take = 6j = q. 
□ 

Proof of Corollary First, we claim that 

max f{x) = max \f(x)\. 

If the degree m of / is odd, this is immediate. If m is even, let us consider the 
polynomial f p for some odd p. Since 

^(• r ) = ai n^*' x ^ wnere a i > o, 

1 iei 



the polynomial f p is also represented as a non-negative linear combination of prod- 
ucts of (q, x), where the cosine of the angle between every pair q, Cj of vectors is 
at least 5. It follows from the proof of Theorem 1.3 above that 



f p dx>0 for any p. 



n-l 



from which we conclude that the maximum value of / and the maximum absolute 
value of / on the sphere S n_1 must coincide. 

Next, as in the proof of Theorem 1.3, we observe that if L C W 1 is a /c-dimensional 
subspace such that for the orthogonal projections c[, . . . ,c' N of ci, . . . , cjv onto L 
we have 

k 

(1 - e)(c-,c') < -(ci,Cj) < (1 - e) _1 (c-,c') for all pairs i,j 

Then 

f / h\ mp / 2 r r 

(1 _ e) m P /2 J fP ^ < / _ \ J fP d[ln < (1 _ e) - mp /2 jf /P ^ 

for all p. In particular, if the degree m of / is even, 

f p dx>0 for all p. 



Therefore, 



s n - 1 nL 



max fix) = max 



The proof now follows from the identities 

l/2p 



/ r \>-/*P 

lim I / f 2p {x) dx ) = max = max f(x), 

p — >+oo \J§n-i J xes™- 1 xes™- 1 

/ /" \ 1/2p 
lim I / f 2p ( x ) dx ) = max \f(x)\ = max f{x), 

p — >+oo \J§n-i n L J xes™~ 1 nL xeS^-inL 

/ f 2p (x) dx = I f 2p dfi n , and 

J sn -i 1 ' 2™ p T(n/2 + mp) J Rn ^ 

[ f 2p (x) dx = -^/^ f f 2p dp k . 

A-inL 2™ p T(k/2 + mp) J L J ^ 



□ 
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3. Examples and an Application 



Some natural examples of sets of vectors q,... ,cjy 6 t" with the property that 
for every the cosine of the angle between q and Cj is at least 5 > are as 

follows. 

(3.1) Examples. 

(3.1.1) Let ci, . . . , cjv G M. n be vectors with positive coordinates such that the 
ratio of the smallest /largest coordinate for each vector q is at least VS. It is easy 
to show that the cosine of the angle between Ci and Cj is at least 5 for each pair 
(hj). 

(3.1.2) Suppose that n = k(k + l)/2 and let us identify R n with the space of 
k x k symmetric matrices with the scalar product (a, b) = trace(a&). Let c\, . . . , cjv 
be positive definite matrices such that the ratio of the smallest /largest eigenvalue 
for each matrix a is at least VS. It is easy to show that the cosine of the angle 
between Ci and Cj is at least 5 for each pair 

Other examples can be obtained by sampling ci, . . . , cn at random from some 
biased distribution in IR n (a distribution with a non-zero expectation). 

Whenever we have a polynomial 

f( x ) = J^J(cj, a;) where ai > 

ic{i,...,n} iei 

|/|=m 

and vectors q as in (3.1.1)-(3.1.2), integration (optimization) of such a polynomial 
over the unit sphere S n_1 reduces to integration (optimization) over a random lower- 
dimensional subspace L. If we want to achieve a (1 — e) m factor of approximation, 
the dimension k of the subspace is only logarithmic in N, so that as long as iV is 
bounded by a polynomial in n, we achieve an exponential reduction in the number 
of variables. 

Finally, we consider the problem of computing (approximating) the hafnian of a 
given positive matrix. This problem is of interests in combinatorics and statistical 
physics and generalizes the problem of computing the permanent, see Section 8.2 
of [Mi78]. Unlike in the case of the permanent, where a polynomial time approx- 
imation algorithm has been recently obtained [J+04], much less is known about 
computing hafnians. 

(3.2) Computing the hafnian of a positive matrix. Let C = (cij) be an 

m x m positive symmetric matrix, where m = 2k is even. Recall (see Definition 
2.1) that the hafnian of C is the polynomial 

haf C = ^ c/, 
i 
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where the sum is taken over all perfect matchings I = ji}, ■ ■ ■ , {ik,jk}j of 

the set {1, . . . , m} and c/ is the product of Cij for {z, j} e /. 

Suppose that C is positive semidefinite. Then C is the Gram matrix of a set of 
vectors, so Cij = (q, Cj) for some vectors ci, . . . , c m e M m and such a representation 
can be computed efficiently (in polynomial time). Using the Wick formula (Lemma 
2.2), we can write 

» m 

haf C = / TT(cj,:r) d/j, m . 

Suppose that for each pair Cj, Cj of vectors the cosine of the angle between Cj and 
Cj is at least 5, which means that c^- > S^/caCjj for every pair Then, by 
Theorem 1.3, to approximate haf C within a factor of (1 — e) m / 2 , we can replace 
the integral by the integral over a random /c-dimensional subspace L C lR m with 
= 0(e _2 5 -2 ln(m + 2)). If e and 5 are fixed in advance, we get a quasi-polynomial 
algorithm of rn°^ nm ^ complexity. 

One can extend the above argument as follows. We observe that haf C does not 
depend at all on the diagonal entries of C, so we are free to change the diagonal 
entries of C to ensure that the above conditions are satisfied. If we put sufficiently 
large numbers on the diagonal of C, we can make sure that C is positive definite, 
so Cij = (q, Cj) for some vectors c±, . . . , c m G M. m . The goal is to make the cosine 
of the angle between every pair Q, Cj of vectors as large as possible. Suppose that 
cu = for all i and let —A be the minimum eigenvalue of C. Then C + XI is a 
positive semidefinite matrix and the cosine of the angle between q and cj is Cij/X. 
Thus as long as the absolute value A of negative eigenvalues of C is sufficiently 
small, we get an efficient algorithm to approximate haf C. 

4. Complex Integration 

Let /, g : M n — >■ R be real n-variate homogeneous polynomials. Let us identify 
W 1 © M n = C n via x + iy = z and let v n be the Gaussian measure on C n with the 
density 

7r _n e - " 2 " , where ||z|| 2 = ||x|| 2 + ||y|| 2 for z = x + iy. 

We recall that z = x — iy is the complex conjugate of z = x + iy. 
Let us define the scalar product on the space of polynomials 

(f,g) = / f{z)g{z) du n 

(although we use the same notation for the standard scalar product on we hope 
no confusion will result since the domains are drastically different). One can easily 
check that the monomials 

x a = x± x . . . x% n for a = (or, . . . , a n ), where > for z=l,...,n. 
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are orthogonal under the scalar product, though not orthonormal: 



(x a ,x^) = 

Therefore, if 



a.\\ . . . a n \ if a = (3 = («i, . . . , a n ) 
ifa^0. 



f =^2 aaX " and # = &aX " 

ceEF ceEG 

are the monomial expansions of / and g, we have 

(/>#) = ^2 acfiaOtil- ■ -a n \. 

aEFnG 

It follows from the integral representation that the scalar product is invariant under 
the action of the orthogonal group: if U is an orthogonal transformation of M. n 
and polynomials /i, g± are defined by fi(x) = f(Ux) and gi(x) = g(Ux), then 
(fi,9i) = (f,g)- 

Various problems of combinatorial counting reduce to computing the scalar prod- 
ucts of two polynomials. 

(4.1) Example. Let a±, . . . , a at and b be some non-negative integer n-vectors. Let 
M be a positive integer. We define 

N 




/w=n E xbi and 9(x)=* b . 



i=l \k=0 



Then the monomial expansion of / contains all monomials x a , where a is a lin- 
ear combination of a±, . . . , ajv with positive integer coefficients not exceeding M. 
Furthermore, if b = . . . , then (/, g) is the number of non-negative integer 
solutions (fci, . . . , k N ), < ki < M, to the equation 

kiai + . . . + kNCiN = b 

times (3\\...(3 n \. The number of such solutions (k±, . . . , /cat) as a function of b 
is often called the vector partition function, cf. [BV97]. Computing the vector 
partition function is generally as hard as counting integer points in a polytope. 

(4.2) Definition. Let us fix a number < 5 < 1 and a positive integer N. We 
say that a pair of homogeneous polynomials /, g : M n — > R of degree m is (5, iV)- 
focused if there exist iV non-zero vectors a±, . . . ,cin G R n and A?" non-zero vectors 
,6tv € M n such that 

• for every pair the cosine of the angle between aj and bj is at least 5; 

• the polynomial / can be written as a non-negative linear combination 



i iei 
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while the polynomial g can be written as a non-negative linear combination 

9( x ) = ^2 Pi 11^'^' 

/ jeJ 

where the sum is taken over subsets I, J C {1, ... , m} of cardinality \I\ = | J| = m 
and aj, (3 j > 0. 

We prove that the value of the scalar product of a well-focused pair of polyno- 
mials can be well-approximated from the scalar product of the restriction of the 
polynomials onto a random lower-dimensional subspace. 

For a /c-dimensional subspace L C M n , let us consider its complexification 
Lc = L © iL C C n . Let Vk be the Gaussian measure in Lc with the density 
n~ k exp{ — ||z|| 2 } for z G Lc- We pick a /c-dimensional subspace L C M n at random 
with respect to the Haar probability measure on the Grassmannian Gfc(M n ) and 
consider the restrictions f^ and ql onto L and the integral 



(/l,#l> = / f{z)g(z) dv k . 

J Lq 

We claim that as long as k ~ log AT, the properly scaled integral over Lc approxi- 
mates the integral over C n within a factor of (1 — e) m . 

(4.3) Theorem. There exists an absolute constant 7 > with the following prop- 
erty. 

For every 5 > 0, for any positive integer N , for any (5, N)-focused pair of poly- 
nomials /, g : M, n — > R of degree m, for any e > and any positive integer 
k > ^e~ 2 5~ 2 ln(iV + 2), the inequality 

1} (f,g)<(i-e)- m (f L ,g L ) 

holds with probability at least 2/3 for a random k-dimensional subspace L C M n . 

The proof is very similar to that of Theorem 1.3. The only difference is that we 
need the complex version of the Wick formula. 

(4.4) Definitions. Let m be a positive integer. A permutation of the set 
{1, . . . , m} is a bijection a : {1, . . . , m} — ► {1, . . . , m}. 

Let C = (cij) be an m x m matrix. The permanent per C of C is defined by the 
formula 

m 

perC = ^ II CiCT W' 

where the sum is taken over all permutations of the set {1, . . . , m}. 

Here is the complex version of the Wick formula. Since the author was unable 
to locate it in the literature, a proof is given here. 
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(4.5) Lemma. Let m be a positive integer and let fi,gi : R n — > R be linear 
functions. Let C = (c^-) be an m x m matrix defined by 



Cij 



= \ fi( z )9j( z ) dvn- 



Then 



P n 

I \~\fi( z )9i(z) du n = perC. 

Jcn i= i 



If fi is defined by fi(x) = (a,i,x) and gj is defined by gi(x) = (bj,x) for some 
cii, bj G R n then = (a,i,bj) . 

Proof. Given vectors a\, . . . , a m and bi, . . . ,b m , let 



P( x ) = Yl(ai,%) and q(x) = 

i=i j=i 

Our goal is to prove that 

{p, q) = per C where c^- = (a^, bj). 

First, we check the identity in the special case when a\ = . . . = a m = e±, the 
first basis vector, and b\ = . . . = b m = b = (Pi, . . . , p n ) is an arbitrary vector. In 
this case, p(x) = x™ and q(x) = (PiXi + . . . + P n x n ) m , so we have (p, q) = P^ml. 
On the other hand, = Pi for all % and j, so perC = miPY 1 as well. 

Next, we check the identity when a±, . . . , a m = a and bi, . . . , b m = b, where 
a and b are arbitrary vectors. Applying scaling, if necessary, we can assume that 
|| a || = 1. Since an orthogonal transformation of R n does not change either (p, q) or 
C, this case reduces to the previous one. 

Now we consider the general case. We observe that both quantities (p, q) and 
perC are multilinear and symmetric in ai, . . . ,a m and multilinear and symmet- 
ric in 61,... ,b m , so we obtain the general case by polarization. For variables 
A = (Ai, . . . , A m ) and [i = (/ii, . . . , (i m ) we introduce vectors a\ = Xiai + . . . + 
A m a m and b^ = ^ibi + . . . + [i m b m . If F(a x , ... ,a m ;bi,... ,b m ) is any polyno- 
mial multilinear and symmetric in ai, . . . ,a m and multilinear and symmetric in 
61, . . . , bm, then (m!) 2 F(ai, . . . , a m ; bi, . . . , b m ) is equal to the coefficient of the 
product Ai • • • A m/ ui • • • \i m in the expansion of F(a\, . . . , ay, 6 M , . . . , 6 M ) as a poly- 
nomial in Ai, . . . , A m , Hi, . . . , [i m . Since if two such polynomials F and G agree on 
all (2m)-tuples (a, . . . , a; b, . . . ,b), they agree everywhere. Letting F = (p, q) and 
G = per C, we complete the proof. □ 



Now the proof of Theorem 4.5 follows the proof of Theorem 1.3. 
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